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Abstract 

We study multi-particle state contributions to the QCD two-point functions of the 
axial-vector and pseudo-scalar quark bilinears in a finite spatial volume. For suf- 
ficiently small quark masses one expects three-meson states with two additional 
pions at rest to have the lowest total energy after the ground state. We calculate 
this three-meson state contribution using chiral perturbation theory. We find it to 
be strongly suppressed and too small to be seen in present-day lattice simulations. 



t Permanent address: Department of Physics and Astronomy San Francisco State University San Francisco, 
CA 94132, USA. 



1 



I. INTRODUCTION 



Hadron spectroscopy is one of the prime applications of lattice QCD simulations. Recent 
results for the light hadron spectrum in 2+1 flavor simulations show an agreement with the 
experimentally measured values to an accuracy of a few percent [1, 2}. The computation 
of the excited-state spectrum is much more complicated, and the numerical results are not 
as satisfactory as those for stable hadrons [3, 4]. Still, steady progress in excited-state 
spectroscopy has been made over the last years (for a recent review see Ref. [5]). 

One feature of unquenched lattice simulations is the presence of multi-particle states in 
the correlation functions measured to obtain the spectrum. In particular, with the up and 
down quark masses getting closer to their physical values, one may expect three-particle 
states with two additional pions at rest to contribute dominantly after the ground state in a 
given channel. Evidence for this expectation has been reported, for example, in Refs. [6, 7]. 
There results for the two-point correlation function C(t) of the pseudo-scalar density were 
analyzed, as a function of the time separation t. The data were well described by an ansatz 
involving two exponentials, 

C{t) = c,e- Mt + Cl e- M '\ (1.1) 

with M' = M + 2m w and the other three parameters determined by the fit. 1 A significant 
dependence of the correlator on the sea-quark mass was observed, supporting the interpre- 
tation of the second term in Eq. (1.1) as a three-particle contribution. 

The question here is whether one really sees a three-particle state and not a genuine 
one-particle excited state. In fact, the results in Refs. [6, 7] are somewhat surprising. On 
theoretical grounds one expects the coefficient c\ of a three-particle state to be suppressed 
by two powers of the (spatial) lattice volume and therefore to be rather small. In principle, 
this volume dependence of the three-particle state can be monitored, providing a strong 
check for the three-particle state hypothesis. However, this would require lattice data at 
various volumes, and these are often not available. 

The main observation in this paper is that the correlation function of the pseudo-scalar 
density and the axial-vector current can be calculated reliably in chiral perturbation theory 
(ChPT) [8-10]. ChPT is expected to give good estimates for the coefficients cq and c\ for 
small pion masses where the exponential suppression of the three-particle state contribution 
in the correlator is less strong. Moreover, to leading order (LO) in the chiral expansion, 
the coefficients c and C\ depend only on the pseudo-scalar masses and decay constants, 
hence they are related. It turns out that, in case of both the axial- vector and pseudo-scalar 
correlation functions, the coefficient c\ of the three-particle state is completely determined 
in terms of c and the Goldstone-meson masses, reducing the number of fit parameters in 
Eq. (1.1) from three to two. 

Our results show that the ratio C\/cq is very small, even smaller than a naive analysis 
would suggest. For typical pion masses and lattice volumes of present lattice simulations it 
is of order 10~ 4 — 10~ 3 . Therefore, the size of the three-pion contribution is smaller than 
typical statistical errors in the lattice data and it plays no role in practice. This seemingly 
negative result can actually be viewed as a positive result, if one is interested in "genuine" 
excited QCD states, and not in the multi-particle state contaminations. 

This paper is organized as follows. After a brief more general discussion of correlators 
in a finite volume (Sec. II), we present our results in Sec. Ill, with subsections III A and 



To be precise, what was analyzed is the effective mass, equal to minus the time derivative of log C(i). 
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Ill B containing technicalities. Section IV confronts our results with lattice data, and Sec. V 
contains our conclusions and some further discussion. 

II. QCD CORRELATORS IN A FINITE SPATIAL VOLUME 

Throughout this article we will consider QCD in a finite spatial box with length L in 
each direction and periodic boundary conditions. The euclidean time extent, however, is 
assumed to be infinite. This choice simplifies our calculations, and is justified as a good 
approximation for the many lattice QCD simulations with T<^2L. 

We are interested in the two-point correlation functions of the (zero component of the) 
axial-vector current and the pseudo-scalar density: 

C a AA (t) = J L3 Sx(A a (x,t)A a (0,0)) , (2.1) 

C a PP (t) = I d 3 x(P a (x,t)P a (0,0)} , (2.2) 

J L A 

where the flavor non-singlet current and density are defined as usual by 

Al = ^, lb T a ^ , P a = i,^T a ^ . (2.3) 

Unless stated otherwise, we consider 2+1 flavor QCD with degenerate up and down quarks 
with mass m and a heavier strange quark with mass m s . Therefore, the flavor index a runs 
from 1 to 8, and the SU(3) group generators T a are chosen equal to A a /2, with A a the usual 
Gell-Mann matrices. 

With our choice for the quark masses the correlators obey isospin symmetry, so it is 
sufficient to study the two cases a = 1 and a = 4. In addition, the correlators (2.1) and 
(2.2) are not independent because the axial-vector current and the pseudo-scalar density are 
related by the partially conserved axial-vector current (PCAC) relations, 

d,Al = 2mP 1 , d,Al = (m + m s )P 4 . (2.4) 

Using these relations the pseudo-scalar correlator can be obtained from the axial-vector one 
by taking two time derivatives. In the following we will thus focus on C% A (t), although much 
of the following discussion applies to C PP (t) as well. 

The integration over the spatial volume in Eqs. (2.1) and (2.2) projects on states with 
zero total momentum. Hence, the dominant contribution comes from the single-particle state 
with the particle being the appropriate pseudo-scalar meson at rest. For the axial- vector 
correlator we obtain 

C a AA , ln (t) = 1*1 |<(M(0)Wp = 0))| 2 . (2.5) 

In order to keep the notation simple we will often refer to all pseudo-scalars as pions and 
denote them by n a , even though for a > 3 they correspond to the kaons and the eta. The 
same applies to the corresponding masses m a , which will denote either the pion, kaon or eta 
mass, depending on the value of the index a. In addition, we will always assume that t > 0, 
so that \t\ = t. 

Note that for Eq. (2.5) we assumed the standard normalization of one-particle states, 

Mq)Mp)) = S ab L 3 2E p - a 8 P4 , (2.6) 
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in a finite volume. Here E$ a = ^p 2 + ml is the energy of the pion 7i a , and 5^ is a 
Kronecker delta. With this convention the state \ir a (p)) and all other one-particle states 
have mass dimension —1. Consequently, the matrix element 

(0\A a o (0)\n a (p = 0)) = m a f a , (2.7) 

defining the decay constant f a in the usual way, and the correlator C a AA l7r have their standard 
mass dimensions of 2 and 3, respectively. 

The axial-vector current excites other states with the same quantum numbers as well. 
The contribution of an excited pion 7i' a has the same form as Eq. (2.5) with the appropriate 
mass m' a ^> m a . However, for sufficiently small pion masses one expects the three-particle 
state with two additional pions to have a smaller energy than this excited state. For this 
three-particle contribution one finds 

CU**(t)=± £ 5 f+ , +l0 i- - \(0\A (0)\nMn(q)n b (k))\ 2 e' E ^ , (2.8) 

p,q,k k,b 

where E tot is the total energy of the state, and we did not show a possible wave-function 
renormalization. For weakly interacting pions E tot equals approximately the sum E^ a + 
E^ b + E^ b of the individual pion energies and we expect the wave-function renormalization 
to be close to one. 

In a finite volume with periodic boundary conditions the momenta p = 27m/ L are quan- 
tized, with n having integer- valued components. If 27r/(m 7r L) is not much smaller than one, 
states with at least one of the pions having a non-zero momentum will be signficantly more 
energetic than those in which all three particles are at rest, and it is sufficient to keep the 
latter contribution only, i.e., 

C a AA , 3n (t) = ^ \{0\M0)K(p = 0)n b (q = 0K(£ = 0))| 2 e -(-+ 2 ^ . (2.9) 

Here the index b is restricted to the values 1, 2, and 3, such that the three-particle state is 
indeed a "pseudo-scalar-plus-two-pion-state." This state is also lighter than a three-particle 
state with two additional kaons. Note that the dimensions in Eq. (2.9) come out right: The 
matrix element on the right hand side is dimensionless, so the left hand side has dimension 
3, as the leading one-particle contribution. 

The three-particle contribution is suppressed by 1/L 6 . This factor is usually taken to 
argue that the three-particle contribution in the correlator is small despite the fact that 
it might be the most important excited-state contribution according to the exponential 
suppression. To make this a little more quantitative consider the ratio of the three-particle 
and the one-particle contributions, assuming on dimensional grounds that the matrix element 
in Eq. (2.9) is of the order m n / f n (for a = 1). We thus estimate 



e~ 2m ^ . (2.10) 



This is a very small number indeed. Typical values for m n L in present day lattice simulations 
are of order 4, while f w L is of order 1.5, leading to a prefactor of the order 10~ 2 or 1CT 3 . In 
the next section more precise results will be obtained within ChPT. 
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III. THE CORRELATORS IN CHPT 



We now turn to a calculation of the ratio (2.10) for various mesonic ground states in 
three-flavor ChPT. We will first set up the framework in Sees. Ill A and III B, after which 
we will present results for the axial- vector correlator in Sec. Ill C and for the pseudo-scalar 
density correlator in Sec. Ill D. In Sec. Ill E we extend our results to include the case of a 
ground-state meson made out of two different quarks both with masses equal to the strange 
quark mass, which is helpful for interpreting some of the results of Refs. [6, 7]. 



A. Basic definitions 

We consider standard finite- volume ChPT for 2+1 flavors in euclidean space-time [10, 11]. 
The leading-order (LO) chiral lagrangian we use reads 

C x = \f tr (fyE^Et) - l 2 fB tr (m+E + Y$M) , (3.1) 
where the field E is the standard SU (3)-valued chiral field containing the pseudo-scalar fields 

S = exp^E^ a ) • (3.2) 

B and / are the familiar low-energy coefficients (LECs) related to the chiral condensate and 
the decay constant in the chiral limit. To LO we have f — f n — fx, and our conventions 
correspond to /„. = 92.2 MeV. The mass matrix M in the mass term of the chiral lagrangian 
is, according to our setup in the previous section, given by 

M = diag(m, m, m s ) . (3.3) 

Expanding the chiral lagrangian in powers of the pion fields we obtain the well-known LO 
pseudo-scalar masses 



m 2 a = < 



2Bm = ml , a = 1,2,3 , 
B(m + m s ) = m\ , a = 4, . . . , 7 , (3.4) 



3 B(m + 2m s ) = 



With the spatial momentum discretized in a finite volume, the (euclidean) pseudo-scalar 
propagators are given by 

G ab (x-y,t-t>) = (n a (x,t)n b (y,t')} = 5>*<*-«>-i- e -^.l*-f| , (3.5) 

with the energy E^ a = y/p 2 + m a . 

The axial-vector current and the pseudo-scalar density corresponding to the definitions 
in Eq. (2.3) are 



A; = l -f tr(T a (Et^E - S9 M Et)) , (3.6) 
p a = l fB tr(T a (E - E f )) , (3.7) 
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Since we will be working to leading order (LO) in this article, we will only need these LO 
expressions. Higher-order terms, involving Gasser-Leutwyler (GL) coefficients [10], will thus 
not be needed. 

Expanding Eqs. (3.6) and (3.7) to first order in the pion fields and substituting into 
Eqs. (2.1) and (2.2) it is simple to obtain 

C a AA {t) = \fm a e- m «\ (3.8) 
C a PP (t) = -(^ e - mat . (3.9) 

As already mentioned, these correlators are not independent, but related by the PCAC 
relations (2.4), which can be used to show that 

Cpp(t) = ~K d 2 tC a AA {t) , (3.10) 

m a 

which is indeed satisfied by Eqs. (3.8) and (3.9). The PCAC relations holds to all orders in 
the chiral expansion and for all n-particle contributions. We will use it in section HID to 
calculate the three-pion contribution in C PP (t). 

A graphical representation of these LO correlation functions is given in Fig. la. The two 
open squares represent the one-pion terms for either the current or the density, which (after 
contraction of the fields) result in the propagator represented by the line. 



B. Multi-meson state contribution to the correlators 

The results (3.8) and (3.9) fall off exponentially with the pion and kaon mass for a = 1 and 
a = 4, respectively. The diagrams in Figs, lb to le lead to contributions with the anticipated 
exponential decay with rates 3m n for a = 1 and rriK + 2m w for a = 4. In these diagrams the 
solid square represents the three-pion terms of the current or density. These terms contain 
three pseudo-scalar fields leading to three propagator lines emanating from them. If two of 
the three propagators correspond to pion propagators the diagram leads to an exponential 
with exponent — (m a + 2m 7r )t. Similarly, the circle represents the insertion of an interaction 
vertex stemming from terms with four pseudo-scalar fields in the chiral lagrangian. The 
vertex includes an integration over the intermediate space-time point where the vertex is 
placed. Also diagram le with two LO currents or densities placed at times t and and two 
inner pion propagators has a contribution with time dependence exp[— (m a + 2m n )t]. 

Note that all diagrams in Figs, lb to le are of the same order in the chiral expansion. In 
order to see this, recall the power counting of finite-volume ChPT in the p-regime, for which 

m quark ~ 0(p 2 ) , 1/L ~ d, ~ 0(p) , 7r a ~ 0(p) . (3.11) 

With these counting rules diagram la counts as order p 4 , since each axial-vector current 
involves one pion field and one derivative. Diagram lb counts as order p 8 , since here each 
current consists of three pseudo-scalar fields and one derivative. For diagrams lc to le 
with one or two vertex insertions, note that each vertex involves four pseudo-scalar fields 
associated with either two derivatives or one power of m quar k. This is of order p 6 . Taking 
into account the integration over space time, which counts as order p -4 , a vertex insertion 
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□ □ 

a) 




e) 

FIG. 1: Feynman diagrams for the axial- vector correlation function. The squares represent the 
axial-vector current at times t and 0, where the open and solid squares denote the one-pion and 
three-pion terms in the current, respectively. The circles represent a vertex insertion at an interme- 
diate space time point; an integration over these points is implicitly assumed. The lines represent 
propagators in position space. At least two of the three lines connecting a solid square or circle 
are pion propagators in our calculation. 

effectively counts as order p 2 . Hence, the diagrams in Figs, lc-e count as order p 8 as well. 
Similar arguments can be made for the pseudo-scalar correlator, and one finds again that 
the diagrams in Figs, lb-e are all of the same order in the chiral expansion. Diagrams 
involving GL coefficients Lj, stemming from higher-order terms in the current and the chiral 
lagrangian [10], are at least of order p 2 higher since they contain either two more partial 
derivatives or one extra power of the quark mass. 

We are interested in the contributions to the correlator that decay with 3m n (a = 1) or 
m K + 2m n (a = 4). This allows us to simplify the calculation by ignoring many terms in the 
expressions for the axial-vector current and the interaction vertices that do not contribute. 
For the correlator C a AA 3n we need to keep only those terms in Aq that involve one 7r a field 
and two light fields, i. e., fields with flavor index 1, 2 or 3. For the interaction vertices we need 
to keep the vertices with at least on n a in order to contract with a field from the current. 
Two more fields need to be light fields in order to end up with two pion propagators in the 
intermediate lines in the diagrams. With these restrictions in mind we can work with the 
following (incomplete!) currents: 

4,nlo = 37(^2^2 - d^l + (2^3))+..., (3.12) 
^Inlo = ^(^TCbd^b- d^n^ +... . (3.13) 
where the dots indicate the terms not needed for our calculation. In the first line one has 
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to add the same terms with the replacement 2 — > 3 for the flavor index, as indicated. In the 
second line a sum over the light flavors is implied (b = 1, 2, 3). 
The relevant interaction vertices needed for diagrams lb-e read 



v 1 = 



--^2 ( " K(d^ 2 ) 2 + 87^x7^2 - 4(^7T 1 ) 2 7T 2 2 + (2^3) 

+ ml(nt + 2vr 2 7r 2 2 + 27r 2 7T 3 2 ) j + . . . , (3.14) 



v 4 = -^772" ( nKd^b) 2 - 27r 4 <9^7r 4 (7r 6 <9 /1 7r 6 ) + (<9 M 7r 4 ) 2 7r 2 



24/2 



+ (m 2 Y + m 2 )7r 2 7r 2 +... , (3.15) 



where again the dots stand for terms not needed for our calculation, and in v 4 a sum over 
b = 1, 2, 3 is implied. 



C. Results for the axial-vector correlator 



The diagrams in Fig. 1 are straightforwardly calculated. Two remarks, however, seem 
to be appropriate. First, recall that we are interested only in the contribution with all 
particles at rest. Hence, although the diagrams lb-e are two-loop diagrams, we do not 
need to perform the sum over internal momenta, cf. Eq. (2.8). Obviously, this simplifies the 
calculation significantly. Second, note that diagram le also contains contributions with a 
time dependence proportional to texp(— m a t). These contributions result in the standard 
renormalization of the pseudo-scalar masses, and can simply be omitted for our purposes. 

For the pion correlator (a = 1) we find 



C\ A (t) = affaire 



1 + 



-2m.i 



512(AL)4(rrvL)2 



(3.16) 



In writing Eq. (3.16), we have made the replacement / — > f n in the three-pion contribution, 
since the difference is of higher order in the chiral counting. Note that the prefactor of 
the three-pion contribution is of the same form as in Eq. (2.10). There is a numerical 
factor 5/512 in the ChPT result which leads to an even larger suppression of the three-pion 
contribution in comparison with our naive estimate (2.10). 
For the kaon correlator (a = 4) we obtain 



c 4 AA {t) 



\f K m K e- mKt 



1 + 



1 



hi 



2A{f K L)\m^Lf L \m K 



-2m 7r t 



(3.17) 



where the function hi(x) is given by 

hi(x) = (1 



x + 



2x 2 



4(1 + a;) 



(3.18) 



As in Eq. (3.16), we have here substituted / — > fx- However, to the order we are working we 
could equally well substitute / — > f n or replace (/i^L) 4 by (/^L) 2 (/ 7r L) 2 in the denominator 
of Eq. (3.17). The difference is of higher order in the chiral expansion than we consider here. 
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FIG. 2: The functions hi(x), % = 1,2 that enter the results for C AA (t) and Cp P (t). 



The function h\(x) is a well-behaved function of the mass ratio m^/rriK, monotonically 
falling from about 0.9 at the physical mass ratio x = 139/493 0.3 to about 1/7 at equal 
masses (see Fig. 2). 

Note that the result (3.17) is not applicable for m n equal or near uik- The reason for this 
is that we omitted the three-kaon contribution to the correlator, because we assumed the 
kaon to be significantly heavier than the pion. Because of this Eq. (3.17) does not reproduce 
Eq. (3.16) for uik = rn^. 

As a technical aside we mention that diagrams lb-d vanish in the calculation of Eq. (3.16). 
The reason is isospin symmetry and the fact that all pions have vanishing spatial momenta. 
The particular form of the current (3.12) implies that diagrams lb-d are proportional to 
energy differences of the pions and thus vanish if all pions have equal spatial momenta. 
Therefore, only diagram le contributes to C AA3n (t). 

We emphasize that the three-pion contribution in Eq. (3.16) is completely determined 
in terms of the leading one-pion contribution: Only f n and appear as non-trivial, un- 
known parameters. In other words, employing Eq. (3.16) in the analysis of actual lattice 
data leads essentially to a two-parameter fit, involving the same unknowns as a standard 
single-exponential fit. The same applies to the kaon correlator once /„- and have been 
determined. 



D. Results for the pseudo-scalar correlator 

The pseudo-scalar correlators are easily obtained using the PCAC relations (2.4). For 
flavor index a = 1 we obtain 



f 2 B 2 

Cpp[t) ~ 2m. 6 



45 

. + 512(Umm„L) 2 6 



-2m 7T t 



(3.19) 



The three-pion contribution is enhanced by a factor 9 compared to Eq. (3.16). This is 
a simple consequence of the two time derivatives in Eq. (3.10) acting on the exponential 
exp(— 3m w t). These result in a factor 9m 2 with the 9 multiplying the factor 5/512. 

The kaon correlator is obtained analogously. In this case the enhancement factor is given 
by the ratio (m^ + 2m 7T ) 2 /m 2 K , so we find 



f 2 R 2 

C 4 PP (t) = -JJ^-e-mxt 



2m 



K 



1 + 



1 



ho 



24(f K L)\m^L) 2 z Vm x 



The function h,2(x) is related to the previously defined hi(x) by 

h 2 (x) = (l + 2x) 2 h 1 (x) , 



(3.20) 



(3.21) 



and sketched in Fig. 2. 

In principle, the enhancement up to a factor of 9 of the pseudo-scalar correlator gives 
some preference to the axial-vector correlator in the determination of the mass and decay 
constant of the pion and kaon. However, the three-particle contribution is very small for 
both correlators so that it does not play a role in practice in any case. 



E. Partially quenched mesons with one or two heavy quarks 

So far we considered 2+1 flavor QCD and ChPT. Useful extensions are the partially 
quenched (PQ) variants in which the masses of some of the valence quarks are chosen to 
be different from any of the sea quark masses [12]. Theoretically, these theories have a 
larger field content due to the presence of additional ghost fields. Although technically more 
involved, the calculation of the previous section can in principle be repeated for the PQ case. 
A full PQ calculation is beyond the scope of this paper, but for this subsection we only need 
to consider two special cases. 

First, we claim that the result (3.17) for the kaon correlator applies without change to 
the Nf = 2 theory with two light quark flavors and a quenched strange quark. Quenching 
of the strange quark is achieved by introducing a ghost quark with mass m s which exactly 
cancels the contribution of the strange quark to the fermion determinant. In PQChPT the 
larger field content appears explicitly because the non-linear E matrix of Eq. (3.2) becomes 
a 4 x 4 matrix. The mass matrix in Eq. (3.3) changes to M — diag(m, m, m s , m s ). 

Now suppose we are interested in the exp[— (uik + 2m v )t] contribution to the axial-vector 
correlation function C\ A . As before we can simplify the calculation of the diagrams lb-e by 
dropping terms in the axial-vector current and interaction vertices that do not contribute 
to the exponential we are interested in. Doing this we end up with the same expressions as 
in Eqs. (3.13) and (3.15): The three-pion term in the axial- vector current needs to contain 
at least one 7r 4 field and two more light pion fields. There are no other terms than those 
already shown in Eq. (3.13) fulfilling this requirement, with or without the ghost fields. The 
same is true for the vertices. Three of the four pion fields in the vertices need to be one 7r 4 
and two light fields, and Eq. (3.15) contains all of those. We conclude that result (3.17) for 
the 2+1 flavor theory applies to the theory in which the strange quark is quenched as well. 

Next, we want to generalize our results to a partially quenched setup with two light flavors 
and two quenched heavy flavors. We will refer to these heavy flavors as strange and charm, 
even though we consider the case with equal quark masses only, taking m c = m s . This 
theory contains a U D S " meson, made of the strange and charm quarks, even though with 
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our choice for the quark masses its mass is not the physical D s meson mass. Nevertheless, 
we are interested in the axial-vector correlation function that projects on this D s meson. 

The full PQ theory would contain two ghost fields compensating the effects of charm 
and strange in the fermion determinant. In the corresponding chiral effective theory the 
field £ becomes a 6 x 6 matrix. However, applying the same arguments as before, one 
can convince oneself that all the additional fields due to the presence of the ghosts do not 
contribute to the terms in the correlator we are interested in. It is therefore sufficient to 
consider standard ChPT with four flavors in order to calculate these terms. In this theory the 
flavor index a in Eq. (3.2) runs from 1 to 15 and the mass matrix (again) assumes the form 
M = diag (m, m, m s , m s ). For the group generators T a we take the obvious generalization 
of the Gell-Mann matrices to the group 577(4). 

Expanding the chiral lagrangian in powers of the pseudo-scalar fields we obtain, in addi- 
tion to the masses in Eq. (3.4), 



B(m + m s ) = m 2 K , a = 9, . . . , 12 , 
2Bm s = m 2 Ds , a = 13, 14 , 
a = 15 . 



_ \B{m + 5m s ) = m~ , 
Note that with our choice for the quark masses we have the relation 



m 2 D = 2m 2 



K 



(3.22) 



(3.23) 



between the LO masses. 

The axial-vector correlation function that projects on the D s meson is the one with flavor 
index a = 13 or 14. To LO we find the results already given in Eqs. (3.8) and (3.9) for the 
correlators. The dominant three-meson contribution is the one with two kaons (defined as a 
meson made out of either the strange or charm quark and one of the light quarks) and one 
pion, as one can easily see from the quark-flow equivalents of Fig. 1. The D s +2n state, on the 
other hand, does not couple to C\\ to the order we are working here. Calculating diagrams 
lb-e and keeping only the contributions with the exponential fall off exp[— (2m K + m n )t] 
we find: 



(t) 



fk m Ds 



-m D t 



1 + 



1 



3(/ 7r L) 4 (m ir L) 2 mn s \m K 



m D , 3 \m K ) 



-(2m K +m 7T -rn Ds )t 



(3.24) 



where the function hs(x) is given by 



hs(x) 



l-x- 



2 - 4x - 2x 2 + x 3 
4(1 + x) 2 



(3.25) 



As before, the pseudo-scalar density correlator is easily obtained from Eq. (3.24) using the 
PCAC relation. Here the enhancement factor of the three-meson contribution is (2m^ + 
m n ) 2 /m 2 Ds , so we find 



ff 



2m 



D c -mp„t 



D s 



1 + 



with h\{x) defined by 



3(f n L) 4 (m n L) 2 m Ds \m K 



h t \ ( 2 + :r ) 2 h i \ 
h 4 {x) = h 3 (x) 



'" :: / /( ( "' ' ) g-(2/»K+/» -m D , U 



(3.26) 



(3.27) 
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FIG. 3: The functions h ifi s(x), i = 3,4 that enter the results for C^ A (t) and Cj? P (t). 



Comparing the results (3.24) and (3.26) with their counterparts (3.17) and (3.20) for the 
kaon correlator there appears to be an enhancement factor 8 for the three-meson contribution 
(1/3 versus 1/24). However, this is not the case when one takes into account all factors. 

Incorporating the factor m 7T /mr, s = Jx 2 /(2 — x 2 ),x = m^/mx, Fig. 3 shows the functions 



hi, eS (x) = ^j—^Hx), i = 3,4, (3.28) 

that enter the results for the correlators. These functions are roughly a factor 4 smaller than 
the ones shown in Fig. 2, effectively reducing the enhancement to a factor of about 2. 

F. A comment on higher-order contributions 

We have assumed that we can take the energy of the three-particle state to be equal to just 
the sum of the individual masses, as would be the case if the three mesons would not interact. 
The pseudo-Goldstone boson character of the pions implies that they interact only weakly 
if all three pions have small (or vanishing) spatial momenta. Interaction corrections thus 
appear only at higher order in the chiral expansion and are suppressed by additional inverse 
powers of the volume. Figure 4 shows an example of a higher order diagram, generated by a 
vertex insertion into diagram lb. It consists of two more propagators and one more vertex 
compared to diagram lb. According to the counting rules of Sec. Ill B, it is suppressed 
by an additional factor 1/L 3 . The time integration leads to term with time dependence 
~ t exp(— 3m n t), which, after exponentiating, corresponds to the expected energy shift of 
order L~ 3 [13]. To the order we are working here such interaction corrections can be ignored. 
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FIG. 4: Higher-order Feynman diagram for the axial- vector correlation function. 



G. An estimate for contributions with non-zero momentum pions 

In Section II we assumed that the three-particle contribution to the correlators is dom- 
inated by the state with all three mesons at rest, cf. Eqs. (2.9) and (2.8). For sufficiently 
large spatial extent L (at fixed pion mass) this is not justified because the "non-zero mo- 
mentum states," i.e., states where some of the pions have non-zero spatial momentum, are 
no longer well separated from the one with all momenta equal to zero. The calculation of 
the contributions due to the non-zero momentum states is beyond the scope of the present 
paper. However, a simple estimate for this contribution will be useful in the next section. 

Let us consider the first subleading contribution in Eq. (2.8) with two of the three pions 
having opposite momenta with magnitude 2n/L. The energy of such a pion can be written 
with the scaling factor 



St) 



4tt 2 

v - \ 1 + 7 TV* , 3.29 



hence E tot = m 7T (2s p + 1) and 8EpE^E^ = 8m^.Sp for the contribution of this state to 
Eq. (2.8). Defining r(t) as the ratio of this contribution to the dominant one in Eq. (2.9) we 
obtain the simple estimate 



It 



r (t) = ^r e -2m.( Sp -i)t _ (3 30) 



s 2 



The 18 in the numerator is the multiplicity of the non-zero momentum state, counting the 
number of different ways two non-zero and opposite momenta can be distributed among the 
three pions. As mentioned before, making L larger at fixed m n the ratio r will eventually 
grow and the non-zero momentum state contribution cannot be ignored. 

In case of the kaon correlator (a = 4) there is a second state where one non-zero momen- 
tum is carried by the kaon. This leads to a somewhat different estimate, but for order of 
magnitude estimates Eq. (3.30) will be sufficient. 



IV. CONFRONTING LATTICE DATA 

Two collaborations [6, 7, 14] provide sufficient information in order to compare our results 
with data. In these references plots are shown for the effective mass associated with pseudo- 
scalar correlation function, with the effective mass defined by 

ma^(t) = ~]ogC% P (t) . (4.1) 

If only the ground state contributed to the correlator, the effective mass would be equal to 
the ground state mass m a , and hence be constant in time. Any deviation from a constant 
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FIG. 5: The effective mass for the 07 data of Ref. [14] (also shown in Fig. 2 of this reference.) 
The solid blue line corresponds to the ChPT predicition for the three-pion state contribution to 
the effective mass. 



value originates from excited and multi-particle states. Our results for the three-particle 
contribution can be written in the generic form 



(t)=m a (l + k a e~^~ m ^) , (4.2) 



with m a and m' a denoting the ground state and three-meson state masses, respectively, and 
a coefficient k a , determined by the results of the previous section. 

Consider first the data reported in Ref. [14]. These data were generated with Nf = 2 
O(a) improved Wilson fermions. One data set, labeled by "07," has a pion mass of 270 MeV 
and a spatial extent L such that m^L = 4. 21. 2 The pion mass is sufficiently light that we 
may expect ChPT to be applicable. The measured pseudo-scalar decay constant is such that 
f n L = 1.58. 3 The combinations m n L and f^L are the only unknowns entering the coefficient 
k\ in the pionic effective mass (a = 1), 

90 

h = 512(^L)4(m 7r L)2 • (43) 
With the particular values given above we obtain 

fci « 1.6 x 10" 3 . (4.4) 



2 The results we need are found in Tables 1 and 18 of Ref. [14]. 

3 Za necessary to obtain this value can be obtained with Eq. (B.l) of Ref. [14]. Note that we converted 
f n L to our convention in which its physical value is 92.2 MeV. 
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This is a tiny number, leading to less than a per mille contribution to the effective mass. 
Figure 5 shows the effective mass for the 07 data set together with the ChPT prediction 
given by Eqs. (4.2) and (4. 4). 4 The latter appears as a straight horizontal line in this plot. 
Apparently, the statistically significant deviation from a constant value in the effective mass 
data cannot be attributed to the three-pion state contribution. 

Let us turn to the data from the CERN group [6, 7]. Figure 2 of Ref. [6] shows the 
effective mass of the pseudo-scalar correlator for two different lattices, labelled D 2 and _D 4 , 
respectively. These lattices were also generated with two flavors of non-perturbatively 0(a)- 
improved Wilson fermions. The measured values for the pseudo-scalar masses and decay 
constants are collected in Tables 9 and 10 of Ref. [7]. 

For the D 2 lattice we focus on the data set with equal sea and valence quark masses. 5 
Here we find f n L = 1.0 and m^L = 5.9. Note that for this data set the pion mass is about 
620 MeV and the applicability of ChPT is questionable. Nevertheless, here we find the value 

fci « 4.4 x 10~ 3 , (4.5) 

which is of the same order as Eq. (4.4). Again, this value is orders of magnitude too small to 
explain the curvature in the left panel of Fig. 2 of Ref. [6]. Strictly speaking, our calculation 
for k\ does not exactly correspond to the lattice data shown in the left panel of Fig. 2, 
because it shows PQ data with the average valence quark mass roughly ten percent smaller 
than the sea quark mass. However, this is a small effect that cannot account for the huge 
discrepency between the data and our ChPT result. 

The right panel of Fig. 2 in Ref. [6] shows the result for the effective mass with two heavy 
valence quarks or, with our nomenclature of Sec. HIE, for m 13cff (t). The data comes from 
the D 4 run with f w L = 0.91 and m^L = 4.1. The D s mass is not given in Ref. [7], but the 
value m Ds = 567 MeV can be read off from Fig. 2. 6 The kaon mass m K = 502 MeV for this 
lattice follows from Eq. (3.23). Using these values in 

1 (2m K + - m Dt )m„ ( m w \ 

we obtain the estimate 

fcia « 2.1 x 10" 2 . (4.7) 

As in the previous cases, this value is too small to account for the curvature in the right 
panel of Fig. 2 of Ref. [6]. 

Result (4.7) for k 13 is larger by a factor 4.7 compared to k\ in (4.5). Most of this 
enhancement stems from the smaller pion mass on the D4 lattice, a smaller part is a partial 
quenching effect. In order to roughly disentangle both effects let us compute k\ for the D4 
lattice: 

ki PS 1.5 x 10~ 2 . (4.8) 

This is already 3.4 times larger than Eq. (4.5). The remaining factor 1.4 compared to 
Eq. (4.7) can be attributed to the valence quark masses being heavier than the sea quark 
masses. 



4 We thank Francesco Virotta and Rainer Sommer for providing us with the lattice data for this plot. 

5 Second row for the D 2 run in tables 9 and 10 of Ref. [7]. Za = 0.75 is given in Ref. [15]. 

6 The figure shows the result for = 0.13590 [16]. 
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Finally, let us estimate the contribution of the non-zero momentum states using 
Eq. (3.30). We choose a sufficiently small time separation t where the exponential fac- 
tor does not completely suppress this contribution and where the data differs significantly 
from the leading ChPT result. 

For the ALPHA data shown in Fig. 5 we choose t/a = 14, leading to r(t) « 1.3. Thus, 
the non-zero momentum state contribution is not negligible. Still, even if we allow for a few 
non-zero momentum states their cumulative contribution is still too small to describe the 
data. The same conclusion can be drawn for the data of the CERN group. If we choose 
t = 0.7 fm we obtain r(t) 1.3 (left panel) and r(t) 0.5 (right panel). 

V. CONCLUDING REMARKS 

Our results make it extremely unlikely that a multi-pion state is seen in pseudo-scalar 
density and axial-vector correlation functions. The overlap of these local operators with 
multi-pion states is tiny, much too small to be seen in current lattice QCD data. This result 
is qualitatively in accord with findings in excited state spectroscopy [3-5], where the data 
showed no indication for the presence of multi-meson states, and extended operators need 
to be employed in order to uncover such states. 

Of course, it is interesting to speculate on how the dependence of the effective mass plots 
of Ref. [6] on the (sea) pion mass might be explained. Possible excited states are the single- 
particle resonance corresponding to the 7r(1300), and, for instance, the p — n or er(500) — n 
excited states. 

The latter two are again multi-particle states, but will be suppressed by only a factor 
linear in the spatial volume, instead of quadratic. Therefore, the factors corresponding to 
k a in Eq. (4.2) will be enhanced compared to the cases we considered in Sec. IV. For the 
pseudo-scalar channel, a p— ir two-particle state has angular momentum L — 1, and therefore 
the p and ir mesons should each have at least one unit of momentum ±2n/L. Using lattice 
data for the vector meson provided in Ref. [7] the two-particle energies are at least about 
2 GeV for both cases considered in Sec. IV, respectively. This is probably too large to fit 
the time dependence of the effective mass plots in Ref. [6]. 

Assuming that the overlap with a er(500) — 7r two-particle state with the local pseudo- 
scalar density is very small, we speculate that the mass m' a in Eq. (4.2) may be a resonance 
for the simulations of Refs. [6, 7}. From the Fig. 2 of Ref. [6] we estimate m' a and k a , finding 
for the left panel m' a ~ 1.7 GeV and k a ps 1.8, and for the right panel m' a ~ 1.5 GeV and 
k a p^ 2. 2. 7 We may fit the dependence on the sea pion mass using the simple form 

m 'a = m 'o + bm l, k a = k + aml. (5.1) 

From such a simple-minded fit we find m' a ~ 1.3 GeV at the physical pion mass, and the 
estimates a ~ —2.0 GeV" 2 and b ps 1.0 GeV" 1 . While we cannot tell from this simple 
exercise whether this is the correct interpretation of the time-dependence of the effective 
masses in Ref. [6], these numbers appear to be quite reasonable, and suggest that the first 
excited state seen in Ref. [6] might correspond to the 7r(1300). Of course, a much more 
detailed investigation, including multiple operators and several spatial volumes, would be 
needed to determine the nature of the time dependence seen in Ref. [6] . 



7 These values for m' a are significantly lower than our estimates for the energies of p — tt states. 
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Finally, we anticipate that conclusions similar to those reported here should also hold 
for other correlaters, such as for instance a nucleon correlator. Chiral symmetry relates 
processes with different numbers of soft pions, and thus we expect that the contribution 
from a state with two extra pions at rest in the nucleon channel can also be estimated using 
the appropriate chiral effective theory. 
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